Line shape of -0(3770) in e+e~ — > DD 
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Interference phenomena observed in the ^(3770) resonance region in the reactions e + e~ — >■ DD 
axe analyzed. To avoid ambiguities in the determination of the 0(3770) resonance parameters, when 
analyzing data between the DD and DDir thresholds, the amplitudes satisfying the elastic unitarity 
requirement should be used. In the absent of information on the P-wave of DD elastic scattering, 
the 0(3770) parameters, determined by fitting the e + e~ — > DD data, can essentially depend on the 
model used for the total contribution of the resonance and background. The selection of the models 
can be toughened by comparing their predictions with the relevant data on the shape of the i/>(3770) 
peak in the non-DD channels e + e~ — »7XcO, ^/V"7> <M, etc - 
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I. INTRODUCTION 



The resonance -0(3770) was investigated in the re- 
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actions e + e — > DD by 
|L MARK-II EL BES 

12, EU, Belle [igj], and KEDR [20K23J Collaborations 
With increasing accuracy of measurements, there ap- 
peared indications on an unusual shape of the 0(3770) 
peak, i.e^, on possible interference phenomena in its re- 
gion [l,[a[lll£ljl3-lll]. Resently, the KEDR Collabo- 
ration noted [2l| - t23l | that the parameters of the ■0(3770) 
resonance become distinctly different from those quoted 
by the Particle Data Group (PDG) in the preceding re- 
views (see, for example, Ref. [23]) if the data analysis 
takes into account the interference between the -0(3770) 
production amplitude and the nonresonant DD produc- 
tion one. In Refs. [22], |23j two very different solutions 
for the interfering resonance and background amplitudes 
were obtained [28j |. These solutions lead to the same 
energy dependence of the cross section and are indistin- 
guishable by the x 2 criterion. Ambiguities of this type 
in the interfering resonances parameters determination 
were found in Ref. (29| . 

CLEO-c has now accumulated about 818 pb -1 (30j 
and BES III about 2.9 fb -1 [3l[ integrated luminosity 
on the 0(3770) peak for open charm physics investiga- 
tions. Therefore, from CLEO-c and BES III, one can also 
expect new data with very high statistics on the shape 
of the -0(3770) resonance. In this regard, we believe it 
is timely to discuss some dangers which are hidden in 
the commonly used schemes for the description of the 
0(3770) peak. 

Section II shows that the most precise current data 
on the e + e _ — > DD reaction cross section in the 0(3770) 
region are hard to describe with a single 0(3770) reso- 
nance. In Sec. Ill, simple models for the isoscalar part of 



the D meson electromagnetic form factor, Fp, which de- 
termines the amplitude e + e~ — >• DD in the 0(3770) res- 
onance region, are constructed. The models take into 
account interference between the resonance and back- 
ground contributions and yield good descriptions of the 
data. The form factor Fp is constructed in such a way as 
to guarantee at least at the model level the elastic uni- 
tarity requirement. Information on the P-wave of DD 
elastic scattering could be a great help in constructing 
the D meson electromagnetic form factor. However, such 
information is not available. Therefore, it is reasonable 
that the -0(3770) resonance parameters, derived from fit- 
ting the e + e~ — > DD data, can essentially depend on the 
model used for the sum contribution of the resonance 
and background. Section IV shows that the selection of 
the models can be significantly toughened by compar- 
ing their predictions with the relevant data on the shape 
of the 0(3770) peak in the non-DD channels, such as 
e + e~ — >-7Xc0j J/ipVi 0*7, e ^ c - The results of our analysis 
are briefly formulated in Sec. V. A comment concern- 
ing the ambiguity of the fitting solutions found in Refs. 
[22l I23I 29] is given in the Appendix. 



II. THE 0(3770) RESONANCE IN e+e" 



DD 



Figure [T] shows the data for the sum of the 
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D + D reaction cross sec- 
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tions in the 0(3770) region, <j(e + e~ —> DD), ob- 
tained by BES H, 11] (68 points in the yfs region 
from 3.645 GeV to 3.872 GeV), CLEO (1 point 
at 3.774 GeV), BABAR 0, II1_(15 points for 

3.73 GeV <yfs < 3.89 GeV), and Belle[19| (9 points for 
3.73 GeV < y/s < 3.89 GeV). Here y/s is the energy in the 
DD center-of-mass system. This is the most detailed and 
accurate current data. It is clear, however, that further 
improvement of the data and matching the results from 
different groups are required. 

Note that the BES Collaboration 0, Q measured, in 
the region up to the DD* threshold (« 3.872 GeV), the 
quantity R(s) = a(e + e~ — > hadrons) / a(e + e~ /i + /i~) 
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Figure 1: The data from BES @, @|, CLEO [Tf|, BABAR 
[13, [3, and Belle QjJ for a(e + e~ — > DD). For clarity, the 
inset shows only the points from CLEO 0, BABAR [rfl.[Isl|. 
and Belle fijl. 



Figure 2: The results of the fit using Eqs. |TJ and J2]). The 
illustration of the 0(3770) resonance shape dependence on the 
parameter r. See the text for details. 



(where tr(e + e 



) =4ira 2 /3s and 



1/137). 



The DD events were not specially identified. The BES 
points shown in Fig. Q] correspond to the cross section 
(Ana 2 /3s)[R(s) - R uds ], where R uds = 2.121 [§] describes 
the background from the light hadron production. This 
cross section gives a good estimate for cr(e + e~ — » Z?_D) in 
the if) (37 70) region, because it is expected that the decay 
width of the cc state "0(3770) (or states) into the non- 
DD modes must be comparable with the decay width 
of the ip(2S) resonance located under the open charm 
production threshold. Hence, the ratio -B(-0(377O) — > 
non-DD)/B(ijj(3770) -t DD) must be small owing to 
the large total width of the ■0(3770). This is confirmed 
by experiment [H, HH . 

The measured DD mass spectrum 0, [lj]-[l9| has the 
following features. First, the right-side of the 0(3770) 
peak turns out to be more steep than its left-side. Sec- 
ond, there is a deep dip near 3.81 GeV in the mass dis- 
tribution (in fact, the cross section dips to zero near this 
point). These features are hard to describe with the help 
of a single 0(3770) resonance contribution. 

In most experimental works, the e + e~ — > DD cross 
section caused by the "0(3770) resonance production was 
described with minor modifications by the following for- 
mula (below for short 0(3770) is also denoted as ip"): 



ay(e + e~ -> DD) 



127riy/ e+e - TjpnjjB ( s ) 

(m%„ sf + (m^rp,( s )y 



(1) 



where my is the mass, T^,,i e + e - the e + e decay width, 
T^„ D £,(s) the DD decay width, and r^'(s) the total 



decay width of the resonance. The energy dependent 
width ^miuq(s) was taken in the form 



(s) = Gl 



P 3 o(s) 



P+(s) 



\l + r 2 p 2 (s) ^ \ + r 2 p\(s)) ' 



(2) 



and 



"D + 



where po(s) 

are the D° and D + momenta, respectively, r is the DD 
interaction radius [33j , and G^,» is the coupling constant 
of the 0" to DD. Because the i\>" — > DD decay is dom- 
inant we put in Eq. (JTJ) T*p(s) =T^„ d£) (s). This 
simplification is not essential for our analysis. 

The dashed and solid curves in Fig. [5] show the fits 
to the data in the region 3.72 GeV < y/s < 3.9 GeV with 
the use of Eqs. (TJ and © at r = and 100 GeV -1 , re- 
spectively. In the inset on this figure the quantity x 2 , 
characterizing the goodness-of-fit, is shown as a function 
of r. As r increases from approximately to 15 GeV -1 
(as 3 fm) the \ 2 value sharply decreases and then, with 
r increasing, remains practically unchanged. Such a be- 
havior of x 2 leaves the parameter r very uncertain. Of 
course, too large values of r hardly have any physical 
means [34j . Owing to the parameter r in F^i, D £)(s), one 
successes in raising the left-side of the -0(3770) peak and 
lowering its right-side. In fact, all existing data require 
such a deformation of the ip" peak. However, as is seen 
from Fig. [TJ a dip near 3.81 GeV cannot be explained by 
varying r. Obtained very unsatisfactory \ 2 values (for 
the dashed and solid curves in Fig. [Ux 2 ~ 413 and 248, 
respectively) are due to both notable differences between 
the data from different groups and existing the dip in 



the DD mass spectrum (for example, for the solid curve 
in Fig. H the points at ^ = 3.8 and 3.81 GeV yield 
X 2 ~ 81). In order to qualitatively improve the data 
description in the tp" resonance region, in particular, to 
explain a dip near 3.81 GcV, it is necessary to take into 
account the interference between the resonant and non- 
resonant DD production. 



III. THE D MESON ELECTROMAGNETIC 
FORM FACTOR 



A. Unitarity requirement 

In constructing the model describing the process 
e + e~ — > DD one must keep in mind that we investigate 
the D meson electromagnetic form factor, the phase of 
which in the elastic region is completely fixed by the uni- 
tarity condition (or the Watson theorem of final-state 
interaction). Experiment clearly indicates that we deal 
with the resonant scattering of D mesons. Really, there 
is the tp" resonance between the DD and DD* thresholds 
(2m D fa 3.739 GeV and m D +m D * pa 3.872 GeV), which 
in a good approximation can be considered as an elastic 
one, because it has not appreciable non-DD decays [28| . 
Usually, such scattering is described as resonance scat- 
tering with an elastic background, see, for example, Ref. 
36], i.e., the corresponding strong amplitude Tj with the 
definite isospin I and spin J (in our case it is the DD 
scattering amplitude T®) is given by [13] 



T? = e i5 ° sin 6° = 



bg _ I 



2i 



2iSb g t* 
c- ■* res j 



(3) 



were 5® = 5b g + S res is the scattering phase, 5b g is 
the elastic background phase (or the phase of poten- 
tial scattering), and S res is the phase of the reso- 
nance amplitude T res (in the simplest parametrization 
T rPS = , f , n ). Then, according to the unitarity con- 

the D meson isoscalar form factor 

si L L 1 



M—E 

dition ImFg, 



r/2' 

has the form in the elastic region 



Figure 3: The graphical representation of the strong DD scat- 
tering amplitude T° and the D meson electromagnetic form 
factor . The vertical dashed lines show that the D and D 
mesons in the loops are on the mass shell. Diagrams corre- 
sponding to the amplitude 75° and the form factor f% show 
the structure of primary mechanisms included in the model 
to describe the tp" resonance region. 



ment, the easiest way to use the field-theory model shown 
in Fig. [3] and write 



where 



T°( S ) 



l-i7?(a) ' 



p o / \ _ fp( s ) 



l-t7?(«) 



l?(s)=v(s)t° 1 (s), 



v(s) = \pl(s)+pl(s)]/VS, 



(6) 
(7) 

(8) 
(9) 



F D =e<G° D = e i ^ +5 ^G° D , 



J D 



(4) 

where G° D is the real function of energy. A similar repre- 
sentation of the amplitude e + e~ — > DD used for the data 
description guarantees the unitarity requirement on the 
model level. The sum of the e + e~ DD reaction cross 
sections is expressed in terms of Fp in the following way 



1 9. 



ip"DD 



6ir m^„ — s ' 



/£(«) = A 7 + 



9ip"~t9iP"DD 



III 



■4," 



(10) 



(11) 



(s) 



8na 2 
3^72 



B. 



Graphically the amplitude Ti(s) and the form factor 
\ f d( s )\ [Po( s ) +P 3 +( s )] ■ ( 5 ) F%(s) defined in Eqs. © and © correspond to the in- 
finite chains of the diagrams in Fig. [3] with the real D 
and D mesons in the intermediate states. The amplitude 
t\ and the form factor /£, defined in Eqs. (JTUJ) and (fTTj) 
specify the structure of primary mechanisms included in 
the model to describe the ip" resonance region. The 
constants A and A 7 effectively take into account back- 
ground (nonresonant in the tp" region) contributions to 
the strong amplitude and form factor, respectively, and 



A simplest model for F^: resonance plus 
background 



To understand how the form factor and strong ampli- 
tude can be constructed satisfying the unitarity require- 



4 



the constants g^n D Q and £y,» 7 describe couplings of the 
ip" to the DD and virtual 7 quantum. The requirement 
of the unitarity condition is fulfilled in the model under 
consideration: the phase of the form factor F D (s) is de- 
fined by the phase of the amplitude T°(s). This phase 
has the dynamical origin. 

The physical content of Eqs. ^ and ([7]) will become 
more clear if they are rewritten in the form of Eqs. ([3]) 
and respectively. As a result, we obtain the follow- 
ing expressions for the background and resonance com- 
ponents of T°(s): 



T bg - 



2iSb g (s) 

2i~ 



1 — iv(s)X ' 



V^ib" Dd( s ) 



Ml„ - s + Re±V (Ml,) - IV {s) 



(12) 



(13) 



where 



ImIV(s) = \T sT 4,"dd{ s ) 



4"dd( s ) 



Gtt 



!/(*), (14) 



ReLV(s) = -A' 



<4>"DD 



67T 



Ka) 2 



ML =ml„-ReU ri (ML), 



9iP'>dd( s ) 
For ^(s) we obtain 



9tp"DD 

|1 - w(s)A| 



(16) 



(17) 



' D[S) " 6 1 |M|„ - s + Ren^(M|„) - LV( S )| ' 

(18) 

where Si(s)=Sb g (s) + S re s(s) [Sb g (s) and S re s(s) are the 
phases of the amplitudes (fT2"j) and (fT3")) . respectively] and 



A 7 (s) 



A, 



\l-iv(s)\\ 



(19) 



Thus Fp(s) incorporates the resonance contribution 
(proportional to <7</>" 7 ) modified (dressed) by the strong 
background [39[ and the proper background contribution 
(proportional to A 7 ) modified by the strong resonance 
and background final-state interactions. The numerator 
in Eq. (|T5| is proportional to the first-degree polyno- 
mial in s, A 7 (m 2 „ — s) + g^n^g^i, £>£>, with real coeffi- 
cients. This ensured that the dip in a(e + e~ —> DD) near 
3.81 GeV can be explained by the zero in Fp(s), caused 
by compensation between the ip" resonance and back- 
ground contributions. Note that the presence of the zero 
in F D (s) is in qualitative agreement with the coupled- 
channel model prediction [40(. 
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Figure 4: The resonance plus background model. The solid 
curve is the result of fitting the data with the use of Eqs. ((5} 
and (|18p . The dashed curve shows the contribution from the 
ip" resonance production (~ gy/' 7 in F%) and the dotted curve 



(15) shows the contribution from the background production ( 



A 7 in F D ) modified by the strong resonance and background 
final-state interactions. 



As is seen from Fig. |4l the constructed model for F^ (s) 
yields the quite reasonable description of the data (here 
X 2 ~ 123, which is much better than the above x 2 va l~ 
ues for the fits shown in Fig. [5]). For the solid curve 
in Fig. |4l the cross section at the maximum (located 
at yfs = ^fs max ~ 

3.773 GeV) 9.13 nb, the full 

width of the peak at its half maximum Th max rs 29.7 
MeV, and the effective electron width of the resonance 
structure T e J + f e _ = s max a max T hmax /(12Tr) « 0.263 keV. 
These characteristics of the observed peak are in close 
agreement with the values of the mass (« 3.773 GeV), 
the total width (« 27.2 MeV), and the electron partial 
width (« 0.262 keV) which are quoted by PDG [H[ as the 
averaged individual characteristics of the ip" resonance. 
However, the peak (in its line shape there is a zero at 
\fs m 3.814 GeV) does not correspond to a solitary reso- 
nance. Therefore, it is reasonable that the model param- 
eters for ip" differ from the effective parameters of the 
visible peak. Let us present the corresponding numbers. 

The curves in Fig. [4] correspond to the follow- 
ing values of the fitted parameters: =3.799 GeV, 



94," dd = ± 19-35, 
GeV" 2 , and A, 



9^- 



= ±0.1483 GeV" 



±25.07 [if 



g^" 1 g^"DD>^ (<o), 

vidual characteristics 



A 7 >0 



, A = -30.35 
(<0), then 
As the indi- 



see Eq. [fig])] 

of the ip" resonance, one can 
take the quantities dressed (renormalized) by the back- 
ground contributions [see Eqs. ([T5 ]) -(fT7 |l ]: Mq,,= 3.784 
GeV, T^p D£> =T^, D o{M^,)/Z^, =37.61 MeV, and 



5 



pren 



- ijj" 



- jZ^n = 0.05181 keV, where Zy> = 1 - 



Ren;„(M2„) = 1.748 and I> e+e - = 4 7 ra 2 3 2 „ 7 /(3M3 / ,). 

The obvious drawback of the considered model is 
the uncertain nature of the background contributions. 
Therefore, the validity of this model is hard to verify 
in other reactions. However, the model can be easily 
improved. It is clear that the main sources of the back- 
ground in the ip" region are the tails from the J/ip, ip(2S), 
-0(4040), ^(4160) and other resonances. The right num- 
ber of resonances can be incorporated in the model by 
adding the corresponding pole terms to expressions (p~0|) 
and (fTTj) for i?(s) and /^(s). In that case, the parameters 
A and A 7 will effectively describe the contributions from 
the residual background, and it is hoped that they will be 
small. The ip(2S) resonance is closest to the ip". Its cou- 
pling to e + e~ is about an order of magnitude larger than 
that of ip" [HI and there are no apparent reasons for the 
suppression of the coupling of the ip(2S) to DD. In the 
next subsection we will consider in detail the model tak- 
ing into account the ip(2S) resonance contribution and, in 
Sec. IV, discuss additional ways of checking this model. 



C. The model for F& with the xp" and -0(25) 
resonances 

The connection of the ip{2S) contribution does not 
change the structure of expressions © and ([7]) for Tj 3 
and Fp. Only the functions t\ and change. Now they 
are given by 



1 9 



ijj{2S)DD 



1 9. 



tf>"DD 



6n m l(2S) 



6tt m\„ 



fUs) = \ 



9iP(2S)j9tP(2S)DD 9ip"l9^»DD 



(20) 



(21) 



>(2S) 



Hereinafter we use the values of ri 
28] and r^ (2S)e+e - = 2.35 keV [28 



= 3.6861 GcV 



From the relation 



I\/,(2S) e + e - = 47TO! ^,(2S) 7 /( 371 V(2S))> W(5 S et 9^(2S) 7 ~ 

±0.7262 GeV 2 . The coupling constant g^,( 2 s)DD i s a lree 
parameters. 

Owing to the common D°D° and D + D~ decay chan- 
nels, the ip" and ip(2S) resonances can transform into 
each other (i.e., mix), for example, ip" — > DD — > ip(2S). 
Therefore, it is very useful to rewrite Eqs. © and (O 
for the amplitude T® and the form factor in terms 
which would reflect this physical aspect of the model, 
and, in particular, introduce the amplitude describing 
the ip" — ip(2S) mixing. 

Let us write the background amplitude in the form 
similar to Eq. (fT2)) 



T b g - 



o2i5 bg {s) _ ^ 
% 



1 — iv(s)X 



(22) 



The amplitude T res , see Eq. ([3]), corresponding to the 
complex of the mixed ip" and ip{2S) resonances, dressed 



by the residual back grou nd, we represent in the following 
symmetric form [35l l4l[ l42j: 

(m|„ - 5)Imn^ (2g) (s) + (m^gg) - s)Imn^(s) 
D 4 ,n{s)D^ {2S) {s) - n2, /i/)(2S) (s) 

(23) 

where D^/^s) and -D,/>(2S)(s) are the inverse propagators 
of ip" and ip(2S), respectively, 

£V(s) =m%„ -s-n^(s), (24) 
A/,(2S) (s) = ml(2S) - s - n V(25) (s) , (25) 



iV'(«) 



fj 



tP"DD 



6tt 1 — iv(s)\ 



v(s) 



(26) 



, , l y^(2S)DD , , , 

T ^^ = ^ l-i v (s)X via) ^ (2?) 

and n^,//^(2S)(s) is the amplitude describing the ip" — 
ip(2S) mixing caused by the ip" — > DD — > ip{2S) transi- 
tions via the real DD intermediate states, 

, v i 9,p"dd9^(2S)dd , > , , 
IVvps) (-) = ^ i - iv{a) x V{S) - (28) 

Note that the phase of T res is defined by that of denom- 
inator in Eq. pS]). 

For the form factor, we get 



D r , («)i>^(2S) (s) - n 2 ( s ) 



, (29) 



where 



K D d(s) = (jn\„ - s)(m 2 (2S) - s)A 7 (s)± 
+9^(2Sh[ D i>" ( s )94>(2S)dd( s ) + n^»v(2S)(s)ff v ,» D £ 1 ( s )] + 

+94>"l[D^(2S){s)9^'DD{ 3 ) + U i>"4>(2S)(s)g4,(2S)DD( s )} 

(30) 

and after cancellations 

K DB (s) = {m%„ - s)(to 2 (25) - s)A 7 (s)+ 
+ (m 2 „ - s)g 4 , (2S)l g^ (2S)DD {s)+ 



+ ( TO ^(25) _ s )5i/'"75'i/i"Dd( S )- 



(31) 



Here Sy,( 2 s)£>fi0) = 9^{2S)dd/\^--w{s)\\\ 9^"dd{ s ) and 
A 7 (s) are given by Eqs. (JT7J) and (Til?)) , respectively. 

The curves in Fig. [5] correspond to the follow- 
ing values of the fitted parameters: m^» = 3.784 
GeV, g v > D D =±13.21, cty» 7 = ± 0.2237 GeV 2 , 
SV(2S)dd = ± 12 -91, A = 26.89 GeV~ 2 , and A 7 = ± 2.456 
[if A 7 >0 (<0), then g^" 1 g^"DD>^ (<0) ancl 
g^(2S) 1 g^(2S)DD < (>0), see Eq. {31])]. Note 
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Figure 5: The model with the ip" and ip(2S) resonances. The 
solid curve is the fit using Eqs. ((5| and (|29[) . The dashed, 
dot-dashed, and dotted curves show the ip" , ip(2S), and back- 
ground production contributions proportional to the coupling 
constants g^"-y, g^{2S)~t> an d A 7 in Eq. ()31[) . respectively. 



that here |A 7 | is about an order of magnitude smaller 
than in the previous case, as qualitatively expected. 
For this fit, \ 2 ~ 125. The form factor has the zero at 
« 3.816 GeV. 

As the individual characteristics of the tp" res- 
onance, one can take again the quantities dressed 
(renormalized) by the background contributions: 
AV =3.789 GeV, Y™ Djj = r^, D5 (M 2 „)/V' =58.03 

r2 



MeV, and r;™ +e - = I> e+e - /Z^„ = 0.2973 keV, 
where = 1 + RelT^, (M|„) = 0.6905 and 

1^,,^- = 47ra 2 g|,„ 7 /(3Af|„). We calculated the 
above parameters with the use of Eqs. (fT2"j) - (fTT| by 
making the substitution 




3.725 3.75 3.775 3.8 3.825 3.85 3.875 3.9 
Vs (GeV) 




3.725 3.75 3.775 3.8 3.825 3.85 3.875 3.9 
Vi" (GeV) 



Figure 6: The predictions of the model with the tp" 
and tp(2S) resonances. (a) The cross sections and (b) 
the phase shifts of DD elastic scattering in the P- 
wave. In (a), the solid, dashed, and dot-dashed curves 
correspond to a(D°D° -» D°D°) = 3tt| sin 5? (s)| 2 /Po(«), 



37r| sin(5 re3 (s)| /po(s), and 37r| sin5(, ff ( s )l /p (s), respectively. 
In particular, the model predicts that a(D°D° 
(5? = 180°) at Vs « 3.846 GeV. 



D U D U 







A ->■ A 



1 9-<p{2S)DD 
67F m l(2S) - S 



(32) 



i.e., in the tp" resonance region, we included in T{, g the 
total background from the amplitude A and the tp(2S) 
contribution, and took into account in T res the tp" con- 
tribution dressed by this total background. For example, 

2 

at the DD threshold, A + ^ — 9 /^ DD a ~ 4 38 QeV" 2 

instead of A f=s —30.35 GeV -2 in the resonance plus back- 
ground model. 

Thus, the fitting of the mass spectrum in e + e~ — > DD 
permits us to determine the resonance and background 
characteristics in specific models. Nevertheless, the infor- 
mation only on the reactions e + e _ — ¥ DD is still lacking 



to give reliable conclusions about the separate compo- 
nents of the reaction amplitude. The performed analysis 
indicates that these components can be very different in 
the different models. On the other hand, it is clear that 
the interference pattern in the tp" region depends on the 
reaction. Therefore, to toughen the selection of the mod- 
els one should compare their predictions with the exper- 
imental data on the mass spectra for several different 
reactions. 

For example, after the fitting of the e + e~ — > DD data 
we all know about DD elastic scattering in the P-wave at 
the model level, see Fig. O Unfortunately, these predic- 
tions are not possible to verify. However, there are many 
other reactions which can be measured experimentally. 
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IV. THE ip" SHAPE IN NON-DD DECAY 
CHANNELS 

Now we apply the last described model to construct 
the mass spectra in the reactions e + e _ — >7XcCb J/^Vi 
(f>ri. In the tp" region, we restrict ourselves to the contri- 
butions only from the ip" and resonances, taking 
into account their couplings to the 7XcO, J/ipVi an d 4>V 
channels in the first order of perturbation theory. 

The cross section for e + e~ — >ab (ab = ~fXcO, J/tpVj 'Pv) 
can be written as 



a ab (s) 



4na 2 k 3 ab (s) 
3s 3 / 2 



\F ab (sW 



(33) 



where fc ab (s) = y/[s - (m a + m b ) 2 ][s - (m a - m b ) 2 ] / (2-^/5) 
and the form factor 



F ab {s) 



Ti-ab(s) 



D^(s)D i , (2S) {s) - Tll„^ {2S) {s) 



(34) 



where 



Tlab{s) = 

0</>(2Sb[A/>" {s)gip(2S)ab + n.0»^(2S)(s)3i/J"ab] + 
±SV>"7 [Dtp(2S)(s)9ip"ab ± n^,//^(2S) (s)flty>(2S)a&] (35) 

and g^,(2S)ab, g^"ab are the coupling constats of the 
ip(2S), ip" to the ab channel. 

Table I presents information about the ip(2S) and ifj" 
resonances in the 7Xrf), J/4>Vi an d 4>f] decay channels 
[H, 03 1") , which we use to construct the corresponding 
mass spectra. The values for g^(2S)ab indicated in the 
table are obtained, up to the sign, from the data on the 
ip(2S) — > ab decay widths by the formula 



r 9^{2S)ab . 3 , 2 



12tt 



(36) 



Note that the available information about the tp" — > 7XcO, 
J/ipVi 4' r l decays are very poor 28]. Data on the mass 
spectra in these channels are still absent. The cross 



Table I: Information about the 4>(2S) and ip" resonances in 
7XcO, J/i/iV, and 4>V decay channels |28l. [4a - t45l ] (errors < 10% 
are not shown). 



ab channel 


7XcO 


J/tpr] 




B(iP(2S) -> ab) 


9.68% 


3.28% 


(2.8±£;°) x 10- 5 


r^(2S)ai>; keV 


29.4 


10.0 


(8.5±i;S) x 10~ 3 


9ip(2S)ab, GeV -1 


±0.25 


±0.22 


±(2.7±° i) x 10- 4 


cr(e + e _ — > ab) at 


72 ±9 


8.6 ±3.9 


3.1 ±0.8 


= 3773 MeV, pb 









sections of the reactions e + e~ — > 7Xc0, JJ^JL 4>V were 
measured by the CLEO Collaboration |43l - l45j at a sin- 
gle point in energy yfs = 3773 MeV (at the supposed 
maximum). Their approximate values are presented in 
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Figure 7: The cross section for e + e~ — > 7XcO- 
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Figure 8: The cross section for e + e — > J/iprj. 
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Table I and Figs. [7H9] by the points with the error 
bars. They allow us to roughly estimate the coupling con- 
stants 5^" 7Xc0 ~ 0.54 GeV -1 , g^"j/^ v ~ 0.053 GeV -1 , 
and gyifa ~ 1-12 x 10 -2 GeV -1 , by using Eqs. ([55])- 
(|35| . to construct the corresponding cross sections. Here, 
as an illustration, we put g^,(2S)ab and g^'ab > 0, and 

5^(2S) 7 /5-i/J"7 < 0- 

The solid curve in Figs. EHS] shows the e + e~ — >7Xc0, 
J/ipr), (jyq reaction cross sections; the dashed and dotted 
curves show the contributions from the ip" and ip(2S) 
resonances proportional to 

9ip"ab 

and 

[94>{2S)jD^ (s) + gV"7 n ^"-0(2S)(s)]ff^(2S)a& 

in Eq. ([55]) . respectively. Note that the cross section for 
e + e~ <frr] is completely dominated by the ip" contribu- 
tion. 

These examples tell us that the mass spectra in the 
ip" region in the non-DD channels can be very diverse. 
Therefore we should expect that the data on such spectra 
will impose severe restrictions on the constructed dynam- 
ical models. 
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Figure 10: The illustration of the ambiguity of the overall 
phase of the e + e~ — > hh reaction amplitude defined in Eq. 
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V. CONCLUSION 

We tried to show that the shape of the ip" resonance 
keep important information about the production mech- 
anism and interference with background. We have con- 
sidered the models satisfying the unitarity requirement 
and obtained good descriptions of the current data on 
the e + e - — > DD reaction cross section, in particular, in 
the model with the mixed ip" an d ip(2S) resonances. 



We have extracted from experiment g 



n/(47T) 



ip(2S)DD 
13. 

Further improvement of the data and matching the 
results from the different groups on the reactions e + e - — > 
DD can result in the crucial progress in understanding 
the complicate mechanism of the ip" resonance formation. 

As we have shown the measurements of the mass spec- 
tra in the ip" region in the non-DD channels, such as 
e + e - — >7Xc0> J/i^Vi 'PVi etc, will also contribute to a 
comprehensive study of the ip" resonance physics and 
the effective selection of theoretical models. 

Additional information about the ip" in the DD mass 
spectra can be extracted, for example, from weak decays 
B — > ip" K and photoproduction reactions at high energies 
jA^ip"A. 
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If the parametrization of the reaction amplitude has no 
clear dynamical justification, it can lead to unexpected 
problems. 

Here we comment on the ambiguity of the interfer- 
ing resonances parameters determination, which has been 
discovered in Ref. [29[ and discussed in Refs. [H, H3l m 
connection with the ip" resonance parameters. 

To illustrate the ambiguity of resonance parameters 
we use a very simple example. Consider the model 
e + e~ — > hh reaction amplitude (where h and h are 
hadrons) involving the resonance and background con- 
tributions |29| 



F(E) 



A x e l 



M-E- iT/2 



B. r 



(Al) 



Here E is the energy in the hh center-of-mass system, M 
is the mass and T the energy independent width of the 
resonance, A x , ip x , and B x are the real parameters. At 
fixed M and T, there are two solutions for A Xl tp x , and 
Br [21: 



(I) A X =A, B X =B, tp z 



(A2) 



(II) A x = V 'A 2 - 2ABT sin tp + B 2 T 2 , B x = B, 

tan (p x — ~ tan ip + BT / (A cost/?), (A3) 

which yield the same cross section as a function of 
energy, a(E) — \F(E)\ 2 , and differ in the magnitude 
and phase of the resonance contribution. For example, 
at M = 3.77 GeV, T = 0.03 GeV, A = 0.045 nb^GeV, 
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(f = 0, and B = 1.5 nb 1/2 , solution (II) gives A x — \[2A 
and <p x = 7r/4. 

For each energy, the two solutions also give the differ- 
ent overall phase, 6 = 5 res + 8b g , of the amplitude F{E). 
For the above numerical example, the phase Sb g corre- 
sponding to solutions (I) and (II) is shown in Fig. [10] 
by the dashed and solid curves, respectively; the phase 



5 res = arctan[ 2 ^ M r _ £ ^ ] is shown by the dotted curve. The 
origin of the rapid change of the phase Sb g (which is ad- 
ditional to S res ) requires a special dynamical explanation 
(for example, the presence of extra intermediate states) , 
for which we do not see at present any reasons. 
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